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We study adiabatic (curvature) and entropy (isocurvature) perturbations produced during a period 
of cosmological inflation that is driven by multiple scalar fields with an arbitrary interaction poten- 
tial. A local rotation in field space is performed to separate out the adiabatic and entropy modes. 
The resulting field equations show explicitly how on large scales entropy perturbations can source 
adiabatic perturbations if the background solution follows a curved trajectory in field space, and 
how adiabatic perturbations cannot source entropy perturbations in the long-wavelength limit. It 
is the effective mass of the entropy field that determines the amplitude of entropy perturbations 
during inflation. We present two applications of the equations. First, we show why one in general 
expects the adiabatic and entropy perturbations to be correlated at the end of inflation, and calcu- 
late the cross-correlation in the context of a double inflation model with two non-interacting flelds. 
Second, we consider two-fleld preheating after inflation, examining conditions under which entropy 
perturbations can alter the large-scale curvature perturbation and showing how our new formalism 
has advantages in numerical stability when the background solution follows a non-trivial trajectory 
in field space. 
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I. INTRODUCTION 

A period of accelerated expansion - inflation - in the 
early universe has become the standard model for the 
origin of structure in the universe. Inhomogeneities in 
the present matter distribution can be traced back to 
quantum fluctuations in the fields driving inflation which 
are stretched beyond the Hubble scale during inflation. 
In the simplest models of inflation driven by a single 
scalar field, these fluctuations produce a primordial adia- 
batic spectrum whose amplitude can be characterized by 
the comoving curvature perturbation 7?,, which remains 
constant on super-Hubble scales until the perturbation 
comes back within the Hubble scale long after inflation 
has ended. 

As soon as one considers more than one scalar field, 
one must also consider the role of non-adiabatic fluctu- 
ations. This can have important consequences, both in 
affecting the evolution of the curvature perturbation (of- 
ten referred to as the 'adiabatic perturbation'), but also 
in the possibility of seeding isocurvature (or 'entropy') 
perturbations after inflation. 

Previous studies have demonstrated that non-adiabatic 
pressure perturbations can alter the curvature perturba- 
tion on super-Hubble scales either during inflation [y||] 
or after ^^- A general formalism to evaluate the cur- 
vature perturbation at the end of inflation in multiple 
fleld models was developed in Ref. Q. In the presence 
of non-adiabatic fluctuations, one must follow the evolu- 
tion of perturbed fields on super-Hubble scales, in partic- 
ular tracking the perturbation in the integrated expan- 
sion ||7|,p|,p|-p^P] , in order to evaluate the large-scale cur- 



vature perturbation at late times 0-|l|,0,|J8l|,|l|,|T| . 

However no similar formalism has been developed 
so far to evaluate the isocurvature perturbation in the 
general case. Instead, isocurvature perturbations have 
been studied in a number of particular models of in- 
flation |lq,0,n2|. These fluctuations typically arise as 
baryon modes (e.g. Il^]) or cold dark matter modes p8|, 
but neutrino isocurvature modes have also been consid- 
ered [19I. Recently, it has been pointed out [EO 21 1 that it 
is rather natural to expect the curvature and isocurvature 
perturbations to be correlated, which yields distinctive 
observational results [E2| , in contrast to the isocurvature 
perturbations usually tested against observations p3| . 

In this paper we will develop a general formalism to 
study the evolution of both curvature and isocurvature 
perturbations in a wide class of multi-fleld inflation mod- 
els by decomposing field perturbations into perturbations 
along the background trajectory in field space (the adi- 
abatic field perturbation), and orthogonal to the back- 
ground trajectory (the entropy field). We allow an ar- 
bitrary interaction potential for the fields, and, although 
we concentrate upon the case of two scalar fields, the gen- 
eral approach can be easily extended to N fields, where 
there will be iV — 1 entropy fields orthogonal to the back- 
ground trajectory. This was done for a specific assisted 
inflation model in Ref. [gij. We will work in the metric 
based approach of Bardeen ||2^ in order to deflne gauge- 
invariant cosmological perturbations, but our formalism 
can also be applied to the study of multiple scalar fields 
in other approaches [p6|-p8[. 

We begin by reviewing the standard results obtained 
in single field models, emphasizing the suppression of 



non-adiabatic fluctuations on large-scales. We then ex- 
tend our analysis to general two-field models, defining an 
adiabatic field and an entropy field, whose fiuctuations, 
though uncorrelated on small scales, may develop cor- 
relations through the subsequent evolution. We present 
two specific models of two-field inflation, one with non- 
interacting fields, the other a model of interacting fields 
which undergo preheating after infiation. 



II. PERTURBATION EQUATIONS FOR 
MULTIPLE SCALAR FIELDS 

We consider N scalar fields with Lagrangian density: 

1 ^ 



C 



-V{^i 



1=1 



and minimal coupling to gravity. In order to study the 
evolution of linear perturbations in the scalar fields, we 
make the standard splitting </J/(t, x) — ^ fiit) + S(pj{t, x). 
The field equations, derived from Eq. (|l|) for the back- 
ground homogeneous fields, are 



(fii + SHipi 



^^. - , 



(2) 



where Vx — dV/dx, and the Hubble rate, H, in a spa- 
tially flat Friedmann-Robertson- Walker (FRW) universe, 
is determined by the Friedman equation: 



H' 



inG 



V{ipi 



1 



^^i' 






(3) 



with a{t) the FRW scale factor. 

Consistent study of the linear field fluctuations Sipi 
requires that we also consider linear scalar perturbations 
of the metric, corresponding to the line elementn 



ds^ = -{l + 2A)dt'^ 



2aB idx^ dt 



+ a^ [(1 - 2^/>)(5y + 2£;,y] dx'dx^ 



(4) 



where we have not at this stage specified any particular 
choice of gauge ||,||,|o|. 

Scalar field perturbations, with comoving wavenumber 
k — 27ra/A for a mode with physical wavelength A, then 
obey the perturbation equations 



5ipj + iH6(p, H — TrSifii 



E 



-2V^,A + ipi 






A + 3TP + —{a^E-aB) 



(5) 



* We follow the notation of Ref. P9| , apart from our use of 
A rather than (j) as the perturbation in the lapse function. 



The metric terms on the right-hand-side, induced by the 
scalar field perturbations, obey the energy and momen- 
tum constraints 



3H 



(^- 



HA]+^ 



iP + Hia^E-aB) =-47rGdp, (6) 
iP + HA = -47rGdq. (7) 



The total energy and momentum perturbations are given 
in terms of the scalar field perturbations by 



(8) 
(9) 






These two equations can be combined to construct a 
gauge-invariant quantity, the comoving density pertur- 
bation lESl 



5p ~ 3HSq 

y^ W>i (Sifj - ipiA ) - (piSifi 



(10) 



which is sometimes used to represent the total matter 
perturbation. 

Because the anisotropic stress vanishes to linear order 
for scalar fields minimally coupled to gravity, we have a 
further constraint on the metric perturbations: 

(o^E-aB) +H (a^E-aB) +i}-A = 0. (11) 

The coupled perturbation equations (|1)-(||) and ( |ll| ) 
are probably most often solved in the zero-shear (or 
longitudinal or conformal Newtonian) gauge, in which 
a'^Ei — aBg ~ pOf . The two remaining metric pertur- 
bation variables which appear in the scalar field pertur- 
bation equation, A^ — <^ and </>£ = ^, are then equal in 
the absence of any anisotropic stress by Eq. ( pi] ) . 

Another useful choice is the spatially fiat gauge, in 
which ipQ = pG,p9]. The scalar field perturbations 
in this gauge are sometimes referred to as the Sasaki or 
Mukhanov variables [ pl| , which have the gauge- invariant 
definition 



Qi = Sipi -t- -fTil^ ■ 
ri 



(12) 



The shear perturbation in the spatially flat gauge is sim- 
ply related to the curvature perturbation, ^, in the zero- 
shear gauge: 

a''EQ-aBQ^a''E-aB + ]-i, = ^^. (13) 

H n 

The energy and momentum constraints, Eqs. (@) and (ra), 
in the spatially flat gauge thus yield 



a 


—AirGem , 








(14) 


'iAQ = 


-AirGSqq , 








(15) 


n in Eq. 


(10), and from 


Eq. 


i) 


we 


have 



The equations of motion, Eq. (ph, rewritten in terms 
of the Sasaki-Mukhanov variables, and using Eqs. (|lj) 
and ( p^ ) to eliminate the metric perturbation terms in 
the spatially flat gauge, become |32|: 



)/ + 3HQi + -^Qi 



E 



y^. 



SttG 



•fiv.1 



.3 



H 



•^i^j 



(16) 



A. Curvature and entropy perturbations 



The comoving curvature perturbation [B3U34] is given 



by 



11 = 11} 



H 



P + P 



-dq 






(17) 



This can also be given in terms of the metric perturba- 
tions in the longitudinal gauge as p9[ 






(18) 



For comparison we give the curvature perturbation on 
uniform-density hypersurfaces. 



p 



(19) 



first introduced by Bardeen, Steinhardt and Turner |3g] 
as a conserved quantity for adiabatic perturbations on 
large scales |3^J§] . It is related to the comoving curvature 
perturbation in Eq. (|l7|) by a gauge transformation 



C = 7^ 



2p 



3{p + p) \aH 



*. 



(20) 



where we have used to the constraint equation (|lj) to 
eliminate the comoving density perturbation, £„. Note 
that TZ and — C thus coincide in the limit fc — > 0. 

Both TZ and — C are commonly used to characterise the 
amplitude of adiabatic perturbations as both remain con- 
stant for purely adiabatic perturbations on sufficiently 
large scales as a direct consequence of local energy- 
momentum conservation [pj, allowing one to relate the 
perturbation spectrum on large scales to quantities at 



the Hubble scale crossing during inflation in the simplest 
inflation models p5 37|. 

A dimensionless definition of the total entropy pertur- 
bation (which is automatically gauge-invariant) is given 

by 

S = h('-P-'-^), (21) 

\p p J 

which can be extended to define a generalised entropy 
perturbation between any two matter quantities x and 

y- 

' Sx Sy^ 



i>xy — ti 



y 



(22) 



The total entropy perturbation in Eq. ( |2l| ) for N scalar 
fields is given by 



S = 



2(V + 3H Y., ^2 ) 5V + 2VY.I M^'^Pi - ^lA) 



3(2y + 3i/Ej'Pj)E/'/i? 



(23) 



where the perturbation in the total potential energy is 
given by W = Yl,i Vvi^Vi- 

The change in TZ on large scales (i.e., neglecting spa- 
tial gradient terms) can be directly related to the non- 
adiabatic part of the pressure perturbation pU 



7^: 



-3iJ-5. 
P 



(24) 



We will thus now consider the evolution of the adiabatic 
and entropy perturbations in both one- and two-field 
models of inflation. 



B. Single field 

Perturbations in a single self-interacting scalar field 
obey the gauge-dependent equation of motion 



Sip + 3HSip - 



K 



-2V^A- 



f 



iptp I 5ip 



A + ii;+'—{a^E-aB) 



a^ 



(25) 



subject to the energy and momentum constraint equa- 
tions given in Eqs. (^"P). 

The scalar field perturbation in the spatially fiat gauge 
has the gauge-invariant definition, Eq. (|2), 



Q^ = Sip 



H 



v-. 



(26) 



For a single field this is directly related to the curvature 
perturbation in the comoving gauge, where the momen- 
tum, 5q = —ipSip, vanishes 



TT TT 

7?. = V + ^Sip = ^Qin ■ 

ip if 



(27) 



It is not obvious that the intrinsic entropy perturbation 
for a single scalar field, obtained from Eq. (23), 



S = 



2V^ 



3ip^{3Hip + 2V^) 



(p lS(p — ipA J — (pStf 



(28) 



should vanish on large scales. Because the scalar field 
obeys a second-order equation of motion, its general so- 
lution contains two arbitrary constants of integration, 
which can describe both adiabatic and entropy perturba- 
tions. However S for a single scalar field is proportional 
to the comoving density perturbation given in Eq. (10), 



and this in turn is related to the metric perturbation, 'J, 
via Eq. dll), so that [51 



S 



K. 



67rG(^2[3iJ(^ + 2V"^ 



r* 



(29) 



In the absence of anisotropic stresses, ^ must be of order 
Aq, by Eq. (pi]), and hence the non-adiabatic pressure 
becomes small on large scales [|[|9 10 . The amplitude of 
the asymptotic solution for the scalar field at late times 
(and hence large scales) during inflation thus determines 
the amplitude of an adiabatic perturbation. 

The change in the comoving curvature perturbation is 
given by 






(30) 



and hence the rate of change of the curvature perturba- 
tion, given by d\nTZ/d\na ^ (k/aH)'^, becomes negligi- 
ble on large scales during single-field infiation. 



C. Two fields 

In this section we will consider two interacting scalar 
fields, (p = ifi and x = y^2- The analysis developed here 
should be straightforward to extend to include additional 
scalar fields, but we do not expect to see any qualitatively 
new features in this case, so for clarity we restrict our 
discussion here to two fields. 

In order to clarify the role of adiabatic and entropy 
perturbations, their evolution and their inter-relation, we 
define new adiabatic and entropy fields by a rotation in 
field space. The "adiabatic field" , a, represents the path 
length along the classical trajectory, such that 



{cos 9)(j) + {sm6)x, 



where 



cosy = 



sm( 



(,2^_^2 



(31) 



(32) 



^2+^2 



This definition, plus the original equations of motion for 
(j) and X, give 



a + 3Ha + K = , 



(33) 



where 



Va = icos0)V^ + {smd)V^ 



(34) 



As illustrated in Fig. m, Sa is the component of the 
two-field perturbation vector along the direction of the 
background fields' evolution. Conversely, fluctuations or- 



X 



Perturbation 




\ 



Background trajectory 



FIG. 1. An illustration of the decomposition of an arbi- 
trary perturbation into an adiabatic {5a) and entropy {5s) 
component. The angle of the tangent to the background tra- 
jectory is denoted by 9. The usual perturbation decomposi- 
tion, along the (j> and x axes, is also shown. 



thogonal to the background classical trajectory represent 
non-adiabatic perturbations, and we deflnc the "entropy 
fleld", s, such that 



Ss = {coa6)Sx — (sin 



(35) 



From this deflnition, it follows that s =constant along 
the classical trajectory, and hence entropy perturbations 
are automatically gauge- invariant ||40|] . Perturbations in 
da, with Ss = 0, describe adiabatic fleld perturbations, 
and this is why we refer to a as the "adiabatic fleld" . 

The total momentum of the two-field system, given by 
Eq. (§), is then 



K. 



'<i>H.t - x^x, 



~a5a , 



(36) 



and the comoving curvature perturbation in Eq. ( fij\ ) is 
given by 



n = ij + H 



(j)S(j) + xSx 



b^ + x' 



i,- 



H 



Sa. 



(37) 



This expression, written in terms of the adiabatic field, 
cr, is identical to that given in Eq. (g7|) for a single field. 
We can also write Eq. (B^ as 



n = (cos^ 6l)7^0 + (sin^ 9)71 



X ' 



(38) 



where we define the comoving curvature perturbation for 
each of the original fields as 



TZi =ip + ^Sipi = ^Qi . 



(39) 



However, even fields with no explicit interaction will in 
general have non-zero intrinsic entropy perturbations on 
large scales in a multi-field system due to their gravita- 
tional interaction, so that TZi for each field is not con- 
served. Although the intrinsic entropy perturbation for 
each field is still of the form given by Eq. (p8|), it is no 
longer constrained by Eq. (|Ij) to vanish as A: ^ 0. This 
is in contrast to the case of non-interacting perfect flu- 
ids, where it is possible to define a constant curvature 
perturbation for each fluid on large scales [|| . 



The comoving matter perturbation in Eq. (10) can be 
written as 



= (T I (5(7 — (tA I — (T(5( 



2VJ 



s, 



(40) 



which acquires an additional term, compared with the 
single-field case, due to the dependence of the potential 
upon s, where 



Vs = (cos6l)T/x-(sin6')y^, 



(41) 



The perturbed kinetic energy of s has no contribution 
to first-order as in the background solution s = 0, by 
definition. 

The total entropy perturbation, Eq. (E3), for the two 
fields can be written as 



S 



(7(5(7 



+ 3Ha^e6s\ . (42) 



3&^{3H& + 2V„) 

X IVcr & ((5(7 - (jA 

Combining Eqs. (p^, ( jicl ) and (^2|), we can write 



S = 



5* VS.. 



^Ss . (43) 



67rG(72[3iJ(7 + 2K] W J 3(72 
Comparing this with the single-field result given in 



Eq. (29), we see that the entropy perturbation on large 



scales is due solely to the relative entropy perturbation 
between the two fields, described by the entropy field ds. 
The change in the comoving curvature perturbation 
given by pUSq] 



n 






r* 



1 



-H 




(44) 



which can be expressed neatly in terms of the new vari- 
ables: 



where 



. iJfc2 2H-^ 



a 



(45) 



(46) 



The new source term on the right-hand-side of this equa- 
tion, compared with the single-field case, Eq. (|3C|), is 
proportional to the relative entropy perturbation be- 
tween the two fields, 8s. Clearly, there can be significant 
changes to 7?. on large scales if the entropy perturbation 
is not suppressed and if the background solution follows a 
curved trajectory, i.e., ^ 7^ 0, in field space |I^. This can 
then produce a change in the comoving curvature on ar- 
bitrarily large scales (i.e., even in the limit /c ^ 0) [p],p9[. 
Equations of motion for the adiabatic and entropy field 
perturbations can be derived from the perturbed scalar 
field equations (0), to give 



a2 



8-6 + 3H(5(7 ^^\'—^y^^-B'^\S(j 



= -2V„A + G 



k 



A + 3ip+—{a''E-aB) 



and 



2{e6sy ~ 2—ess , 

a 



5s + 3H5s + [ — + Vss-e^ ]Ss 



(47) 



-2- a{6a-&A) 
a 



aSa 



where 

Vss- 



(sin^ 9)Vx^ + (sin26')V0x + (cos^ ^^4,4, ■. 
(sin^ 9)V^^ - (sin26')V0x + (cos^ ^Wxx ■ 



(48) 



(49) 
(50) 



When 9 — 0, the adiabatic and entropy perturbations de- 
coupled. The equation of motion for 5a then reduces to 
that for a single scalar field in a perturbed FRW space- 
time, as given in Eq. (p5|), while the equation for 5s is that 



^If we employ the slow-roll approximation for the back- 
ground fields, (j) ~ —V^/3H and x — ~^x/3^, we obtain 
9 ~ 0. This reflects the fact that the rate of change of 9 is 
slow - instantaneously it moves in an approximately straight 
line in field space. But the integrated change in 9 cannot 
in general be neglected. Even working within the slow-roll 
approximation, fields do not in general follow a straight line 
trajectory in field space. 



for a scalar field perturbation in an unperturbed FRW 
spacetime. 

The only source term on the right-hand-side in Eq. ( p^ ) 
for the entropy perturbation conies from the intrinsic 
entropy perturbation in the cr-field. From Eqs. (IIj) 
and (E^) we have 



&{Sa — a A) — aSa — 2&95s 



AttGo? 



*, 



(51) 



and hence we can rewrite the evolution equation for the 
entropy perturbation as 



6s + "iHSs 



n2 



Vss + W^ I Ss 



& 27rGa^ 



*. (52) 



Note that this evolution equation is automatically gauge- 
invariant and holds in any gauge. On large scales the 
inhomogeneous source term becomes negligible, and we 
have a homogeneous second-order equation of motion for 
the entropy perturbation, decoupled from the adiabatic 
field and metric perturbations. If the initial entropy per- 
turbation is zero on large scales, it will remain so. 

By contrast, we cannot neglect the metric back- 
reaction for the adiabatic field fluctuations, or the source 
terms due to the entropy perturbations. Working in the 
spatially flat gauge, defining 



= dan — Sa 



H 



^, 



and using 



Aq = 47rG— Q^ 



(53) 



(54) 



we can rewrite the equation of motion for the adiabatic 
field perturbation as 



Q, + 'iHQa 




Vaa 



8ttG 



H 



2{e6sy -2\ Y^ + ^\e5s. 

' cr H j 



Qa 



(55) 



When = 0, this reduces to the single-field equation of 
motion, but for a curved trajectory in field space, the 
entropy perturbation acts as an additional source term 
in the equation of motion for the adiabatic field pertur- 
bation, even on large scales. 

In order for small-scale quantum fluctuations to pro- 
duce large-scale (super-Hubble) perturbations during in- 
flation, a field must be "light" (i.e., overdamped). The 
effective mass for the entropy field in Eq. (p2|) is /i^ = 
Vss + 3^^. For /ig > |-ff^, the fluctuations remain in the 
vacuum state and fluctuations on large scales are strongly 
suppressed. The existence of large-scale entropy pertur- 
bations therefore requires 



^l = Vss + 39^ < 2-ff' • 



(56) 



III. APPLICATION TO ENTROPY/ADIABATIC 
CORRELATIONS FROM INFLATION 



Equations (52) and (|55| ) are the key equations which 
govern the evolution of the adiabatic and entropy per- 
turbations in a two field system. Together with con- 
straint equations (^fj) and (jsj) for the metric pertur- 
bations, they form a closed set of equations. They allow 
one to follow the effect on the adiabatic curvature per- 
turbation due to the presence of entropy perturbations, 
absent in the single field model. This in turn will allow 
us to study the resulting correlations between the spec- 
tra of adiabatic and entropy perturbations produced on 
large-scales due to quantum fluctuations of the fields on 
small-scales during inflation. 

A useful approximation commonly made when study- 
ing field perturbations during inflation, is to split the 
evolution of a given mode into a sub-Hubble regime 
{k > aH), in which the Hubble expansion is neglected, 
and a super-Hubble regime (fc < aH), in which gradient 
terms are dropped. 

If we assume that both flelds (p and x ^tre light (i.e., 
overdamped) during inflation, then we can take the fleld 
fluctuations to be in their Minkowski vacuum state on 
sub-Hubble scales. This gives their amplitudes at Hubble 
crossing (fc — aH) as 



Qi 



Hh 



J.H 



■.ei{k), 



(57) 



where I = (j),x, Hk is the Hubble parameter when the 
mode crosses the Hubble radius (i.e., Hk = k/a), and 
60 and e^ are independent Gaussian random variables 
satisfying 

(e/(fc)) = , {ei{k)e*j{k')) = 6ij 5{k - k') , (58) 

with the angled brackets denoting ensemble averages. It 
follows from our deflnitions of the entropy and adiabatic 
perturbations in Eqs. ( pl| ) and (35) that their distribu- 
tions at Hubble crossing have the same form: 



Qo 



Hi, 



%H 



=ecr(fc), 6s\ 



Hi, 



k=aH 



■■es{k), 



(59) 



where Ca and Cs are Gaussian random variables obeying 
the same relations given in Eq. (|58|), with I, J = a, s. 

Super-Hubble modes are assumed to obey the equa- 
tions of motion given in Eqs. ( |55| ) and (p2[), which we 
will write schematically as 



0''iQ^) = S^Ss), 



(60) 
(61) 



where 0'^{Qa) and 0^{5s) are obtained by setting /c = 
on the left-hand side of Eqs. (pa) and (p2h respectively. 



and S" (Ss) is given by the right-hand side of Eq. ( p5| ) . As 
remarked before, there is no source term for 6s appearing 
on the right-hand side of Eq. ( |5^ ) once we neglect gradi- 
ent terms. The general super-Hubble solution can thus 
be written as 



Q, = A+/+(i) + A_/_(t) + P(i), 
Ss = B+g+{t) + B^g^it) , 



(62) 
(63) 



where the real functions f± and g± are the grow- 
ing/decaying modes of the homogeneous equations, 
0''if±) = and d''{g±) = 0, and P(t) is a particular 
integral of the full inhomogcneous equation (|60|). Note 
that the growing-mode solution /+ ex &/H. 

Henceforth we shall consider only slow-roll inflation 
where the evolution can be approximated by first-order 
equations [dropping Ss and Qo- in Eqs. ( p2| ) and (p5|)], so 
that we haven 



Q, ~ Af{t) + Pit) , 
,5s ~ Bgit) . 



(64) 
(65) 



We can, without loss of generality, take f = 1 ~ g and 
P = when k = aH , so that the amplitudes of the 
growing modes at Hubble-crossing are given by Eqs. (|5£ 
as 



A{k) 



H,, 



r{k), B{k)^ 



v/2fc3 



es{k). (66) 



From Eq. (BOp , we see that the amplitude of the partic- 
ular integral P{t) at later times will be correlated with 
the amplitude of the entropy perturbation, _B, and we 
can write Pit) = BP{t), where P{t) is a real function 
independent of the random variables 6,7, es. 
In order to quantify the correlation, we define 



{x{k)y*{k')) 



2^ 

fc3 



C,ySik-k') 



(67) 



The adiabatic and entropy power spectra are given by 

2 

f2 



'Pq. = c, 



Q^Qa 



Vss = Css 



27r 



\t 



IP" 



\9 



while the dimensionless cross-correlation is given by 
Cq^Ss _, gP 



(68) 
(69) 

(70) 



■'■We note that in non-slow-roll scenarios the decaying 
modes may not be negligible on super-Hubble scales, which 
could affect the correlations between adiabatic and entropy 
perturbations. 



Note that the adiabatic power spectrum at late times is 
always enhanced if it is coupled to entropy perturbations 
[i.e., P(i) 7^ 0, in Eq. (|64|)], as the entropy field fluctu- 
ations at Hubble-crossing provide an uncorrelated extra 
source. 

As an illustration, we consider the correlations in the 
adiabatic and entropy perturbations at the start of the 
radiation era, produced after double inflation, as studied 
in Ref. [M. The double- inflation potential for two non- 
interacting but massive scalar fields is: 



y-\< 



1 2 2 



(71) 



Following ||ll||, it is possible to parametrise the back- 
ground scalar field trajectory in polar coordinates when 
both fields are slow-rolling: 



X- 



N 
2^ 



sm a , 



N 
2^ 



■ cos a . 



(72) 



where N — — ln(a/aond) is the number of e-folds until 
the end of infiation. The background trajectory can then 
be expressed as: 



7V~iVo 



(sin a 



|2/(fl= 



(cosa)2«V(fl^-i) 



(73) 



where R — vfi^jvfi^. The scalar field position angle, a, 
can be related to the scalar field velocity angle, 0, which 
we used to define the adiabatic and entropy perturba- 
tions: 



tan( 



iV 



iHa V 27rG 



: tan a . 



(74) 



The scalar field x is assumed to decay into cold dark 
matter while the scalar field </> decays into radiation. 
The entropy/isocurvature at the start of the radiation- 
dominated era is described by 



S, 



ad 



Spc 
Pc 



3 (5p-y 
477 



(75) 



In Rcf . [ pO[ , it is shown how the super-Hubble perturba- 
tions in the radiation era can be determined in terms of 
the perturbations during the inflationary era. The fluc- 
tuations in both (f> and x fields can contribute to both 
the adiabatic and entropy perturbations. The adiabatic 
component comes directly from the comoving curvature 
perturbation, TZ, at the end of inflation, and is given by 

T^rad — ~v47rGW— g-Ffc [{s\nak)e^{k) + (cosafc)e0(fc)] . 

(76) 

The isocurvature perturbation at the start of the 
radiation-dominated era is related to the entropy pertur- 
bation between the two flelds at the end of inflation Il2] 

'Sx_S_ 
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(77) 



which yields 



5'rad ^ -V4:7rG\ ^Hk [R*secak + cosecafe] es(fc) , 



(78) 



and 
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e<T(fc) - 
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1 
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The entropy perturbation during the radiation era only 
depends on the entropy perturbation at Hubble-crossing 
during the inflationary era, while the adiabatic perturba- 
tion during the radiation era depends on both the adia- 
batic and entropy perturbations at Hubble-crossing. This 
is consistent with equations ( p2^ and y^), showing that 
the entropy perturbation sources the adiabatic perturba- 
tion on super-Hubble scales, but not vice versa. 

As both equations ( |78| ) and (^ depend on the random 
variable Cs , the adiabatic and entropy perturbations will 
be correlated, and we find 



C 



"^rad 5'rad 
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(80) 



This correlation is investigated fully in [gO| in terms of the 
usual scalar field perturbation variables. An interesting 
point that can easily be seen from Eq. (JT^) is that 7?.rad 



will depend only on eo- if R 
will be no correlation if R 
Eq. (73), a will be constant for R 



/rric/, — 1. Thus, there 



1. 



As can be seen from 
— 1 and thus so will 
6] a straight-line background trajectory will be obtained 
for i? = 1. This is consistent with Eq. (pTJ), where it 
can be seen that the entropy component only sources the 
adiabatic component on large scales if ^ 7^ 0. 



IV. APPLICATION TO PREHEATING AFTER 
INFLATION 



Preheating typically amplifies long-wavelength modes 
preferentially. As discussed in |^,^,^, amplification 
of super-Hubble modes does not lead to a violation of 
causality, due to the super-Hubble coherence of the infla- 
ton oscillations set up by the prior inflationary phase. If 
TZ is amplified on super-Hubble scales, this will alter the 
resulting imprint on the anisotropics of the cosmic mi- 
crowave background (CMB), and break the simple link 
between CMB observations and inflationary models. 

We consider first the case where the inflaton is mas- 
sive (to y^ 0) and neglect its self-interaction (A = 0). 
The traditional resonance parameter for the strength of 
preheating at the end of inflation is 



9M 
Am? 



(82) 



In the massive case, where modes move through the res- 
onance bands of the Mathieu chart, and for inflation 
at high energies where the expansion of the universe is 
very vigorous, q needs to be much larger than one if the 
parametric resonance is to be efficient p3[. It is possi- 
ble to have large q even for small coupling, g^ ^ 1, as 
TO ^ 00 "^ A/pi. We can write the effective mass of the 
X during inflation as 



TO^ 

7J2 



3g Mj 



(83) 



where 0o is the initial value of at the beginning of 
preheating. It then follows from Eq. (p3) that x must be 
heavy during inflation for this simple potential if efficient 
preheating is to be obtained. 

Any change in the curvature perturbation TZ on very 
large scales must be due to the presence of non-adiabatic 
perturbations. In p3,E5|, it was shown how, if m? ^ m?> 
during inflation with A = = 5, then the x field and 
hence any non-adiabatic perturbations on large scales are 
exponentially suppressed during infiation, and no change 
to TZ occurs before backreaction ends the resonance. 

However, when 5^0, the x field will have a nonzero 
vacuum expectation value (vev) during infiation even 
along the valley of the potential. In the slow-roll limit 
for 0, this vev is determined by V^^^ =0, which gives 



In this section we use the entropy/adiabatic decom- 
position of the perturbation equations to investigate the 
dynamics of super-Hubble perturbations during a period 
of preheating at the end of inflation. We consider three 
models, encompassed by the general effective potential 



2 4 



+ ^5'0V+5'0'X- (81) 



The essence of preheating lies in the parametric amplifl- 
cation of fleld perturbations due to the time-dependence 



of their effective mass, e.g. 
simplest cases, the inflaton ( 
of inflation. 



K 



g 4> . In the 



< simply oscillates at the end 



X~ — 2^ 
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(84) 



The g coupling has the effect of rotating the valley of the 
potential - which the attractor trajectory approximately 
follows - from x = 0, through an angle 



9^ 



(85) 



where, to ensure that the chaotic inflation scenario is not 
drastically altered, we assume Efl] 



-<1 
9 



(86) 



The effect of g is to change the attractor for both x 
and 5x during inflation, since the x ^"^^ ^X equations of 
motion gain inhomogeneous driving terms proportional 
to g^4)^ . This does not necessarily imply that TZ will be 
amplified by preheating at the end of inflation as purely 
adiabatic perturbations along the slow-roll attractor now 
have a component along x 9'S well as <j). In order to de- 
termine whether or not the evolution of the comoving 
curvature perturbation, TZ, on super-Hubble scales is af- 
fected, we need to follow the evolution of the entropy 
field perturbationPI, defined by Eq. (|35|), which gives 



5x + 



~9' 



(87) 



In the limit g/ g — > we recover 5s -^ 5x- Crucially, the 
evolution equation (|5^) for the entropy perturbation has 
no inhomogeneous terms in the long-wavelength (fc — > 0) 
limit, even for g ^ 0, and entropy perturbations will only 
be non-negligible on super-Hubble scales if the entropy 
field is light during inflation. 

In the slow-roll limit and on large scales, the evolu- 
tion equation (p3) for the entropy perturbation has the 
approximate solution E7| 



5s 



where 



-3/2 



k 



n 

H^ 



(89) 



and the effective mass of the entropy field, fig is defined in 
Eq. (pq). The power spectrum of entropy perturbations 
is 



P.oci/^l- 



3-2Ro(i/) 



(90) 



The real part of ly vanishes for fJ.^/H^ > 9/4, leaving 
a steep k^ blue spectrum, which is exponentially sup- 
pressed with time. 

Using Eqs. (|^), (|l|), (|85|), and (||), one finds that 
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(91) 



fij./H^ has a local minimum for g = 0. Thus the addi- 
tional g term in Eq. (p|) serves to increase the entropy 



mass relative to the Hubble parameter, and so does not 
avoid the suppression of the entropy perturbation. The 
g term therefore does not significantly alter the spectral 
index of the spectrum of entropy perturbations, which 
remains steep if g 3> 1. The strongly blue spectrum im- 
plies that non-linear backreaction is dominated by small- 
scale modes, which go nonlinear long before the cosmo- 
logical modes, implying that resonance ends before TZ 
changes 



3 44| 



We have also integrated the field equations numerically 
to avoid relying on any slow-roll-type approximations. To 
numerically evaluate the entropy perturbation, one could 
simulate the original perturbation variables dcj) and Sx, 
using Eq. (|^), and then work out Ss algebraically via 
Eq. (35). However, this approach is prone to numerical 
instability when the entropy perturbation is suppressed. 
To illustrate this, we take g = 8x lO^'^g and g = 3.8 x 10^ 
After about 60 e-folds of inflation, one can see analyti- 
cally that Ss r^ 10"^". Numerically, Sxcosd ~ 5(j) sin 9 ~ 
10^* during infiation. So in order to obtain a high enough 
accuracy to model the suppression of Ss, we require that 
5x cos 6 and S4> sin 9 have to be simulated to a relative 
accuracy of ~ 10~^/10~^'^ = lO"'^^. This means approx- 
imately 32 significant figures are needed, which is beyond 
the capability of standard numerical ordinary differential 
equation integration routines. 

If instead we use the new adiabatic and entropy field 
perturbations and integrate Eqs. ( p2[ ) and (pa), then this 
numerical instability does not occur, since one no longer 
needs to find the difference between two nearly equal 
quantities. Simulation results using these equations are 
compared with the results using the old field perturba- 
tion equations (m in Fig. 0. The simulations show that 
the gro wth in TZ is driven by Ss, in concordance with 
Eq. ([45|). As can be seen, the numerical result using the 
field perturbation equations fails to track the exponential 
decay of the entropy during infiation and thus underes- 
timates the delay in the growth of TZ. 

In practice, we find a similar instability if we try to 
construct the gauge-invariant metric perturbation, ^, re- 
quired in Eq. (p3) in terms of the constraint Eq. (|5l|). 
This includes the intrinsic entropy perturbation in the a 
field, which does become small at late times/large scales, 
but results from the diminishing difference between finite 
terms. It is more stable numerically to follow the value 
of ^ at late times using the evolution equation 



* 



H 



* = AttG&Q, , 



(92) 



^ From Eq. mS) we see that dSs must be non-zero to change 
TZ on large scales. Because 9 ~ 0, from Eq. (pq), the entropy 
remains decoupled from the adiabatic perturbation during 
slow-roll inflation in this model. But at the end of inflation, 
during preheating, 6 j^ 0. 



which can be obtained from the definition of ^ given in 
Eq. (nj) and the metric constraint equations (M) and (pd]). 
Note that the adiabatic/entropy decomposition be- 
comes ill-defined if cr = 0, i.e. both fields stop rolling, 
and this can cause numerical instability during preheat- 
ing if the trajectory is confined to a narrow valley. This 
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FIG. 2. Numerical simulations of the entropy and comov- 
ing curvature perturbations during inflation and preheating, 
with A = 0, g = 2 X 10"^, g = 8x lO'^g and m = IQ-^Mpi. 
The 'new' prefix indicates that the field perturbations were 
evaluated by numerically integrating Eqs. (p2), and (p5[), 
while the 'old' prefix indicates that the perturbations were 
evaluated by integrating the original field equations (H). We 
have not included any higher-order corrections such as backre- 
action from small-scale perturbations which would shut down 
the resonant amplification of Ss at some point. 



can occur, for instance, when 5 = and only the (j) field 
oscillates. The original field perturbations Scf) and Sx re- 
main well-defined, although the comoving curvature per- 
turbation TZ, defined in Eq. (37) becomes singular when 
(T = H. This does not happen for the simulation re- 
sults shown in Fig. with g ^ where the fields oscillate 
in a two-dimensional potential well. 

The massive inflaton potential {m ^ 0) safeguards the 
conservation of 7?, by a bootstrap effect: if preheating 
is strong, g ^ 1, then the entropy perturbation is heavy 
during inflation; on the other hand, if the entropy is light 
during inflation, then 9 < 1 and preheating is very weak. 
This is not altered by a rotation of the trajectory in field 
space [g ^ 0) as can be most quickly seen by noting, 
from Eqs. ^ and @, that 



v., 



v.. = y^. 



V. 



XX 



(93) 



Thus if the x field is very massive (V^x ^ H^)^ '^e must 
have Va-a + Ks ^ H^ ■ For slow-roll inflation we require 
Vera < H'^ and hence Vss > H'^- 

This situation does not hold if the entropy mass dur- 
ing inflation is not linked to the entropy mass during 
preheating pq |, or in massless {m — 0) self-interacting 
(A ^ 0) inflation models [|8|,|8l||]. This latter class of 



models are almost conformally invariant, allowing ana- 
lytical results from Floquet theory to be applied. The 
Floquct index, fi^, which determines the rate of expo- 
nential growth, can reach its maximum as k/aH — -> 0, 
when (?^/A = 2n^ for integer n, thereby implying max- 
imum growth for the longest- wavelength perturbations. 
Assuming slow-roll inflation driven by y « A0^/4, we see 
from Eq. (|9|) that Vaa + Vss > V^^ = g^cj)"^ and thus that 
the entropy fleld is massive {Vss > 9fl"^/4) whenever 
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A M? 



(94) 



PI 



However, we can have resonance at large scales for ji ~ 1 
and g^/X — 2, when the entropy field need not be 
heavy during inflation and no exponential suppression 
takes place, so that the subsequent growth of TZ is ex- 
plosive iQ . The growth of TZ occurs before backreaction 
can shut off the resonant growth of the entropy perturba- 
tions Ss |[48| , P8| , pO| , H ■ Although the region of parameter 
space around g^/X = 2 is thus ruled out, the same does 
not hold for g^/X ^ 1, since the entropy field is then 
heavy during inflation and Ss is again suppressed. 



V. CONCLUSIONS 

We have introduced a new formalism in which to fol- 
low the evolution of adiabatic and entropy perturbations 
during inflation with multiple scalar fields. We decom- 
pose arbitrary field perturbations into a component par- 
allel to the background solution in field space, termed 
the adiabatic perturbation, and a component orthogo- 
nal to the trajectory, termed the entropy perturbation. 
We have rederived the field equations in terms of these 
rotated fields in Eqs. (B3) and (pq). These show that 
the adiabatic perturbation on large scales can be driven 
by the entropy perturbation, while the entropy pertur- 
bation itself obeys a homogeneous second-order equation 
on super-Hubble scales. There can only be significant 
change in the large-scale comoving curvature perturba- 
tion if there is a non- negligible entropy perturbation, and 
if the background trajectory in field space is curved. 

Our formalism can be applied to evaluate the correla- 
tion between the adiabatic and entropy perturbations at 
the end of inflation. As an example we considered the 
example of two non-interacting fields in double infiation, 
calculating the cross-correlation between the adiabatic 
and entropy perturbations. 

The effect of preheating on the large-scale curvature 
perturbation can also be addressed within our formal- 
ism. The mass of the entropy field during inflation is a 
crucial quantity. If the entropy field is heavy, then any 
fluctuations on large scales are suppressed to negligible 
values at the beginning of preheating. This squeezing 
of the entropy perturbation is most accurately modelled 
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numerically using our evolution equation for the entropy 
perturbation. If it is estimated from the usual field equa- 

tions it may eontain larpe numerical errors when there 
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